A two-channel R-matrix analysis of magnetic field induced Feshbach resonances 
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A Feshbach resonance arises in cold atom scattering due to the complex interplay between several 
coupled channels. However, the essential physics of the resonance may be encapsulated in a simplified 
model consisting of just two coupled channels. In this paper we describe in detail how such an 
effective Feshbach model can be constructed from knowledge of a few key parameters, characterizing 
the atomic Born-Oppenheimer potentials and the low energy scattering near the resonance. These 
parameters may be obtained either from experiment or full coupled channels calculations. Using R- 
matrix theory we analyze the bound state spectrum and the scattering properties of the two-channel 
model, and find it to be in good agreement with exact calculations. 
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I. INTRODUCTION 



A Feshbach resonance can occur in collision events, 
when a closed channel bound state of a specific electronic 
level of an atom is coupled to the scattering continuum 
of an open channel of another electronic level of the atom 
. If the energy of the discrete closed channel state is 
very close to that of the scattering particles, the cross 
section is resonantly enhanced, due to the temporary 
capture of the particles in the quasi-bound state. Fes- 
hbach resonances would be of limited importance in cold 
atom scattering if not for the ability to control the en- 
ergy of the bare resonance state, and hence the resonance 
position, through the application of a static magnetic 
field 0, Q • In this context the open and closed channels, 
in the separated atom limit, correspond to pairs of al- 
kali atoms in different combinations of eigenstates of the 
hyperfine and Zeeman Hamiltonian with total angular 
momentum projection quantum numbers rriF,i and 777,77.2- 
The entrance channel is in an s-wave orbital state (rela- 
tive angular momentum quantum number 1 = 0), while 
the closed channels may be in any partial wave, with mi 
indicating the projection of the relative orbital angular 
momentum. All channels have the same projection of the 
total angular momentum along the magnetic field axis, 
772 p = mp,i + 777 772 + ™z , and since their magnetic mo- 
ments are different, they experience a relative Zeeman 
shift in an applied magnetic field. Hence a closed chan- 
nel bound state may be shifted in energy relative to the 



open channel threshold, bringing it into resonance. The 
coupling between the entrance channel and s-wave closed 
channels arises from the central part of the molecular po- 
tential, which only depends on the magnitude of the total 
electron spin. Therefore, the interatomic potential has 
short-range off-diagonal terms proportional to the differ- 
ence between the triplet and singlet potentials. Closed 
channels of higher angular momentum (I > 0) can be 
coupled to the entrance channel through the spin-dipole 
inteaction Q. 

The Feshbach coupling gives rise to a dispersive vari- 
ation with magnetic field of the scattering length, well 
described by [3( 



a(B) = a bg 
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This is associated with a bound molecular state of the 
coupled channels problem crossing threshold at B = Bq. 
The width of the resonance is characterized by AB, and 
abg is the limiting value of the scattering length for mag- 
netic fields far removed from Bq. In ultracold atomic 
gases the s-wave scattering length gives a complete char- 
acterization of the interactions, and Q therefore implies 
an unmitigated control of the sign and strength of the ef- 
fective interparticle potential. This has facilitated the ex- 
perimental investigation of the crossover between loosely 
bound Cooper pairs and a Bose-Einstein condensate of 
true dimers in a trapped Fermi gas 0,0, HSU E3 

The coupled channels description of atomic scattering 
may be solved in full generality once the adiabatic Born- 
Oppenheimer (BO) curves for the singlet and triplet po- 
tentials are known, and gives accurate predictions for the 
location of the experimentally observed Feshbach res- 
onances. However, there is decided usefulness in the 
availability of simplified model descriptions, which can 



2 



be included effortlessly in many-body 0, 0, 0, Q 
and three-body |l5j treatments, or are more amenable 
to analysis of time-dependent experiments 0, HfiL . 
Fortunately, in most atoms Feshbach resonances occur 
at well separated magnetic fields. Hence, a natural start- 
ing point for a simplified description is to consider an 
isolated resonance. 

In this paper, we demonstrate that it is possible to re- 
duce the full multi-channel, low-energy scattering prob- 
lem of two alkali atoms to an effective two-channel prob- 
lem for a single Feshbach resonance 0, ^| . The model 
consists of a single closed channel containing a bound 
state, the bare resonance state, which interacts with the 
scattering continuum in the open channel. In this sense it 
is equivalent to Fano's theory of resonances arising from 
a discrete state embedded in the continuum , but the 
bare resonance state in the two-channel model is a con- 
struct, which may represent the interactions between sev- 
eral closed channels in the full coupled channels problem. 
This model is capable of capturing the relevant parts of 
the full molecular physics of the resonance, provided it 
is constructed to comply with a small set of low energy 
parameters. Section 1TT1 gives an explicit prescription for 
constructing the two-channel pseudopotentials. 

The R-matrix formalism, which is presented in Sec- 
tion II I II is an elegant and efficient method of calculat- 
ing both bound state and scattering properties. It ex- 
presses the logarithmic derivative of the wavefunction at 
some cutoff radius, beyond which the asymptotic solu- 
tions of the Schrodinger equation are known, in terms of 
the solutions of an eigenvalue problem. In Section llVl we 
discuss briefly the specific numerical implementation of 
the method we have used in this paper. We argue that 
a Discrete Variable Representation of the Hamiltonian 
combined with a grid based on a Finite Element partin- 
ioning of space defines a numerical procedure, which is 
efficient in the evaluation of matrix elements and capa- 
ble of handling the two widely separated length scales 
associated with the molecular and low energy scattering 
physics. 

Section [V] gives the results of our calculations for a 
particular Feshbach resonance. In addition, we give the 
parameters defining our two-channel model for several 
other resonances of interest. In all cases we have ensured 
that the two-channel model gives an accurate represen- 
tation of the molecular physics, by comparing with full 
coupled channels calculations. We conclude our findings 
in Section lVT1 where we also contrast our pseudopotential 
approach with that of other authors. 



II. THE TWO-CHANNEL MODEL OF 
FESHBACH RESONANCES 

In the case of an isolated Feshbach resonance it is pos- 
sible to reduce the full coupled channels problem to an ef- 
fective two-channel model with one open channel (1) and 
one closed channel (2) with a bound state, which interacts 



with the open channel scattering continuum to produce 
the resonance. Such a model is capable of capturing the 
essential physics of the full multi-channel scattering cal- 
culation over a wide range of energies. In this section we 
set up the two-channel model, following the procedure 
outlined in Refs. [Til Il7|. and discuss how a correspon- 
dence may be established between its parametrization 
and low energy scattering observables, available either 
through experiments or more exact calculations. 

The coupled two-channel Hamiltonian describing the 
relative motion of the particles is chosen to be 



H 



2B — 



2fX dr' 2 

W(r) 



W(r) 



(2) 



where fi is the reduced mass. We take the zero of en- 
ergy to correspond to the energy of the dissociated atoms 
in channel I at all magnetic field values. Correspond- 
ingly, the only magnetic field dependent component of 
the problem is the energy offset of the second channel. 
In the absence of any coupling between the channels the 
eigenstates of the two-channel problem is given by the 
diabatic states (<fii in , 0), and (0, 02, n), where (j>i in are the 
eigenstates of the Hamiltonian for channel i. The cou- 
pling mixes these states, and the resulting dressed (or 
adiabatic) states will in general have components in both 
channels. Our intention is to solve the coupled channels 
eigenvalue problem 



(H 



2B 



EI) 



0, 



(3) 



for the components of the wavefunction in each of the 
two channels, both for the discrete part of the spectrum 
and for the scattering states. Here I is the two by two 
identity matrix. 

The BO potentials are dominated by the van der Waals 
interactions at large internuclear separation, vanishing as 
r -6 , asr-> oo. At short distances there is a steep poten- 
tial barrier. Any potential we choose with these charac- 
teristics will give an accurate representation of the molec- 
ular physics at low energies. For simplicity we choose a 
Lennard-Jones potential: 



Vi(r) =-4e 



(?)'-(?) 



12 



(4) 



As in Refs. 0, 0] we take the potential in both chan- 
nels to have the same form, but with the closed channel 
potential shifted upwards in energy, such that its thresh- 
old is at -E t h + A/X.B, where denotes the difference 
in magnetic moment between the separated atoms and 
the bare resonance state. This simplifies the numerical 
analysis, and does not seem to have a noticeable effect 
on the results. Hence, we have 



V 2 (r;B) = V 1 (r) + E th + A f iB. 



(5) 



Finally, we represent the coupling between the two chan- 
nels by an off-diagonal potential W(r). Again, the spe- 
cific form of this potential is arbitrary, but it should be 
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a short-range function, since the inter-channel coupling 
is only significant at small distances. Specifically, we use 
an exponential form 



W(r) = W Q e- r/a ", 



(6) 



where ao = 0.529177 x 10~ 10 m is the Bohr radius. We 
thus end up with a five parameter model for the Fes- 
hbach resonance: e and a characterize the depth and 
shape of the diagonal potentials, respectively, while the 
off-diagonal coupling is of strength Wo- Finally, the off- 
set between the threshold energies of the two channels for 
a given magnetic field strength is determined by A/i and 
E t h- These five parameters are determined by matching 
the low energy scattering properties of the two-channel 
model with those of full coupled channels calculations 
or experimental observables. The model potentials are 
illustrated in Fig. the case of 40 K. 
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FIG. 1: Diagonal and off-diagonal potential curves for the 
two-channel model representation of the 40 K Feshbach reso- 
nance at B = 202.15 G (1 G=0.0001 T). The open poten- 
tial, Vi(r) is plotted in blue, while the closed channel poten- 
tial at zero field, Vi(r, 0), is represented by the red curve. 
The diagonal potentials are indistinguishable on this scale. 
The interchannel coupling, W(r), is plotted in green. The 
inset shows the long-range behavior of the potentials, and 
the threshold energy is evident as the offset between the 
asymptotic values of the open and closed channel potentials. 
1 E h = 4.359744 x 10" 18 J. 



To make this connection explicit we make use of ana- 
lytical results for low energy scattering in inverse power 
law potentials derived by Flambaum et al. 0]. For the 
open channel, which determines the far from resonance, 
or background, scattering properties, we require that our 
model potential V\ (r) supports iV s bound states and that 
its scattering length equals the background value, ab g , 
which appears in (Q. We only require that N s is of the 
order of the number of bound states in the actual po- 
tential, since when N s 3> 1, the near-threshold proper- 
ties are essentially independent of its actual magnitude. 
Given the atomic mass and the Cq coefficient, which is 



well known from the asymptotic scattering properties, we 
may unambiguously construct a model potential meeting 
these constraints. Both N s and abg are related to the 
WKB phase $ evaluated at zero energy: 



a bg = a [1 
N, 



tan($ - tt/8)] , 
5" 
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1, 



(7) 
(8) 



where [ ] is the integer part, and the mean scattering 
length (or characteristic scale length of the van der Waals 
potential) a — y/ifj,CQT(j)/T( j)h is governed solely by 
the asymptotic behavior of the potential. The second 
equation is necessary to completely determine the semi- 
classical phase, since the first equation can only be solved 
for <J> modulo 7r. Following Flambaum et al. we can then 
find the relationship between the C\i and C§ coefficients 



Cl9 = 



/ 2^r(i)c 6 5/6 ' 
ior(f)<M 



(9) 



Knowing both the Cq and C\2 coefficients of our model 
potential, the range and depth parameters are uniquely 
determined: 



C\2 

4a 6 ' 



i/(> 



(10) 



(11) 



The bare resonance state is a particular unit normal- 
ized bound state of closed channel potential: 



n 2 d 2 , ; 
W + y2(r;0) 



res (r) = £: rcs (O)0 rcs (r). (12) 



We arbitrarily choose (\> ves {r) to be the second to last 
bound state of the closed channel potential, such that 
E Tes (0) = E^+Eth- Here E^\ indicates the binding en- 
ergy of the bare resonance state relative to the threshold 
in the closed channel. This is calculated by diagonalizing 
— (h 2 /2fj,)(d 2 /dr 2 ) + Vi(r). A calculation of the scattering 
properties of the two-channel model 0, reveals that 
the magnetic field dependent scattering length (Q) can be 
related to the strength of the model coupling potential 
via 



AB=lim-™ 
fe^o 2kabgA.fi 



(13) 



where E = h 2 k 2 /2fi. For E > the decay width of the 
quasi-bound molecule giving rise to the resonance is given 
by HEl 

2 



T(E)=2ir\(4> tM \W\i>E)\ A 



(14) 



Here <Pe{t) is the energy normalized outgoing wave scat- 
tering state at energy E in the presence of the open chan- 
nel potential alone 



w +Fl(r) 



4>E{r) = E(j} E (r). 



(15) 
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Its asymptotic form determines the background contri- 
bution to the scattering phase shift, 6b g (E) 



2fi 
irh 2 k 



sin (kr + 5 hg (E)) . 



(16) 



The background scattering length follows from ab g = 
— ta,nSb g /k in the limit of E — > 0. We specify be- 
low how the matrix element of the coupling operator 
in (|13[l is calculated in the R-matrix formalism. Since 
the scattering length a(B) is zero at the magnetic field 
B = Bq + AB, the width of the resonance can be either 
measured or obtained from multi-channel calculations. 
Hence the strength of the interchannel coupling func- 
tion may calculated using l|13|) , assuming the difference in 
magnetic moment of the two channels A[i is known. We 
will explain below how this quantity may be obtained. 

The measured or calculated position of the Feshbach 
resonance, B , is shifted away from the field value B les , 
where the bare resonance state crosses zero energy, due to 
the coupling between the channels . An accurate analytic 
estimate of this shift is given by multi-channel quantum 
defect theory [2^ 



Bo — B ICS — AB 



x(l — x) 

i + (i - xy 



dbg/o 



(17) 



This provides a link to the threshold energy, E t h- Since 
the variation of the bare resonance state with magnetic 
field is given by E YCS (B) — E Tes (0) + AfxB , it follows from 
E Tes (B ieB ) = that 



E, 



th 



(A^S r . 



E, 



(18) 



Finally, we address the derivative of the bare resonance 
state energy with respect to the magnetic field. From the 
Fano theory of resonant scattering we know that a 
bound state embedded in the continuum gives rise to a 
resonant contribution to the phase shift 



5 Tes (E,B) = tan 



T(E)/2 



E - E Ies (B) + A(E) • 



(19) 



The shift of the resonance position away from the energy 
of the bare resonance state is Q, ^| 



A(E) = V 



dE' T(E') 
2tt E-E' 



(20) 



Importantly, this shift is independent of the magnetic 
field. Hence, we expect that for magnetic fields suffi- 
ciently far from B that the resonance energy is unaf- 
fected by threshold effects, the position of the scattering 
resonance may be fitted with a linear function of the mag- 
netic field, with a slope of A/z. This defines a procedure 
for determining the last parameter of our two-channel 
model. 

As we shall demonstrate below, this procedure for con- 
structing an effective two-channel model yields predic- 
tions, which accurately reproduce the results of a more 
rigorous calculation for both the continuum and the 
bound state spectra. 



SOLUTION USING THE R-MATRIX 
METHOD 



III. 



The R-matrix method |20l |21| has been demonstrated 
over the past thirty-five years to be a general and quite 
powerful ab-initio approach to a wide class of atomic and 
molecular collision problems. The essential idea is to di- 
vide space into two or possibly more physical regions. In 
each of these regions the time independent Schrodinger 
equation may be solved using techniques designed to be 
optimal to describe the important physical properties of 
that region. The solutions and their derivatives are then 
matched at the region interfaces. A simple example to 
illustrate this idea is electron scattering from atomic tar- 
gets. In this case, when the electron is close to the other 
atomic electrons, there are strong electrostatic, exchange 
and correlation effects. Outside this internal region, the 
electron being scattered is subject to purely long-range 
electrostatic forces. It then becomes possible to expand 
the wavefunction in the internal region in a discrete set of 
eigenstates using the full power of standard bound state 
approaches. The internal wavefunction is then matched 
to known or easily computed solutions of the long-range 
potentials on the surface of the sphere bounding the two 
regions. 

For the two channel problem considered in this paper, 
using the definitions in Section[Hj we define the R-matrix 
eigenstates, {ip c ,n( r ) ( c = 1) 2)}, inside a sphere of radius 
of r = a, 



(H 2B +LI- E n I) 



where L, the Bloch operator |23j is, 



(21) 



(22) 



The function of the Bloch operator is to ensure that any 
non-zero contributions to matrix elements arising from 
the kinetic energy operator on the surface of the R-matrix 
sphere are exactly cancelled by contributions from the 
Bloch operator. This ensures that the operator H + L 
is Hermitian and that l|21|l has only real eigenvalues. 
We use a notational convention where the bound states 
(E < 0) have n < 0, while the discrete continuum states 
are labeled with n — 0, 1, 2, . . ., their number limited by 
the finite size of the basis. The R-matrix boundary a is 
chosen large enough that the influence of the potentials 
at r = a is negligible, such that the R-matrix eigenstates 
match unto free particle solutions. These cigenfunctions 
are then used to expand the solution to the two channel 
problem at an arbitrary energy, E, as 



4>2,E{r) 



V'i,™( r ) 



(23) 



The coefficients, determined by matching the so- 

lution and its first derivative to appropriate asymptotic 
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forms at r = a. The boundary a is chosen sufficiently 
large that when both channels are closed, corresponding 
to true vibrational bound states of the coupled system 
with E < 0, the wavefunction has decayed to zero at the 
boundary. This ensures that the R-matrix eigenstates 
reproduce the bound state vibrational spectrum of the 
molecules precisely. For the scattering states the match- 
ing conditions give, 



4>2,E{a) 



where 



ft 2 



^-22 



j>c,n(a)lpc> ,n( a ) 
E n — E 



^2,E( a ) 



c=l,2, 



(24) 



(25) 



is the R-matrix for our problem, and tp' c E (a) is the 
derivative of the channel wavefunction at the boundary. 
It is convenient to define the matrix 



,_i /2 _ /_2/WlA/fcT o 



(26) 



which determines the energy normalization of the scat- 
tering states. Here 



fci = 



2(j,E 



k 2 = 



2[i(E - E th - A[iB) 



(27) 



In the case where channel 1 is open and channel 2 is 
closed, the asymptotic scattering state for s-wave colli- 
sions can be written as 



"01, -EM 




[sin(fcir) + K cos(fcir) 
Aexp(— nr)/r 



(28) 

where k = ik 2 , A is an irrelevant normalization constant, 
and the phase shift is given by t&n S(E) = K. Upon use 
of the R-matrix matching condition, the phase shift in 
the open channel is found to be 



tan<5(£T) 



kilZn cos(fcia) — sin(fcia) 
cos(fcia) + kiTZn sin(/cia) ' 



(29) 



where we have defined an effective open channel R-matrix 
element, TZu — TZh — \TZi 2 \ 2 k/(1+TZ 22 k). In the absence 
of coupling, W = 0, we have TZu = TZu and Vi e{ t ) = 
Mr)- 

When both channels are open, we must allow for the 
possibility that the particles may enter the collision re- 
gion in one channel and exit in the other. This is de- 
scribed by a 2 x 2 K-matrix, K, and the asymptotic scat- 
tering state becomes in a general matrix notation 



-1/2 



[f(r) + g(r)K]. 



(30) 



The off-diagonal elements of ipE( r ) describe the cases 
where the incoming channel and outgoing channel differ. 
The matrices /y (r) = sin(kir)Sij andgjj (r) = cos(kir)5ij 



give the regular and irregular solutions in the asymp- 
totic region [24]]. The R-matrix matching at r = a then 
gives a linear system of equations for the multi-channel 
K-matrix: 

[k- 1/2 g(a)-72.k" 1/2 g'(a)]K = [^k- 1/2 f'(a)-k- 1/2 f (a)]. 

( 31 ) 

The S-matrix can be found from 



I + iK 
I -iK' 



(32) 



The multi-channel S-matrix satisfies the unitarity condi- 
tion£JSy| 2 = lforall j. 

A great advantage of the R-matrix method can be 
seen from (|25|) . The explicit energy dependence enables 
the R-matrix to be calculated quite efficiently with en- 
ergy. In the presence of resonances, the phase shift varies 
rapidly with energy but the corresponding R-matrix is 
constructed easily from energy independent parts. In 
addition, the R-matrix poles are often very good approx- 
imations to the positions of the resonances. The widths 
may be extracted in a variety of ways, the simplest being 
computing the phase shift as a function of energy and fit- 
ting the result to a Breit-Wigner form plus background 
term. This makes R-matrix theory ideal for the study of 
magnetic field induced Feshbach resonances. 

The original derivation of the R-matrix method by 
Wigner and Eisenbud [2{| used R-matrix states satisfy- 
ing a fixed, zero derivative boundary condition at the 
R-matrix radius. While this is perfectly acceptable for- 
mally, it is a very slowly convergent process numerically 
as a consequence of the fact that the scattering wavefunc- 
tion at energy E satisfies a different boundary condition 
at the R-matrix radius. Thus, the use of fixed boundary 
condition states, can only be correct in the limit of a com- 
plete expansion basis. The most appropriate way to look 
at the problem from a variational viewpoint is that we 
are finding the Rayleigh-Ritz solution to the self-adjoint 
operator H+L. By chosing a basis which does not satisfy 
a fixed derivative condition at the surface, convergence of 
the solution is quickly achieved with an incomplete basis. 
It is possible to devise an approximate correction proce- 
dure (2(| to account for the infinity of neglected levels 
using a fixed boundary condition basis. In the approach 
above, the use of the Bloch operator eliminates the need 
for specifying a specific derivative condition on the basis. 
In fact, basis functions with arbitrary derivatives may be 
used as an expansion set. This leads to a very rapid con- 
vergence of the S-matrix as a function of the basis set 
size. The higher lying levels have little physical signif- 
icance individually, but their spectral sum is effectively 
converging to the R-matrix. In addition, it can be shown 
that truncating the very high lying R-matrix levels has 
a negligible effect on the low-energy scattering proper- 
ties. This is in contrast to some other approaches [27j 
where the elimination of a single state gives completely 
erroneous results. 
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IV. NUMERICAL APPROACH 

In the R-matrix method as outlined in the last section, 
the critical element is the computation of the R-matrix 
eigenvalues and eigenvectors. The Bloch operator for- 
malism enables the use of powerful variational methods, 
which do not require the underlying basis set to satisfy 
any specific boundary condition. This produces a very 
rapid convergence of the results with respect to basis set 
size. The choice of primitive basis in which to expand the 
R-matrix eigenstates is, of course, arbitrary. In practice 
many choices are eliminated on the basis of computa- 
tional convenience or efficiency. For the problem under 
consideration in this paper, there is the additional con- 
straint that the basis accurately describe both the two- 
body vibrational bound states as well as the very low 
energy continuum states. For the calculations presented 
here this requires a highly accurate description of the 
wavefunction to internuclear distances of approximately 
1000a . The Finite Element Discrete Variable (FEDVR) 
basis, described in other publications |28ll2{^ . has proven 
ideal for this purpose. While a detailed description of the 
method must be relegated to the references, a few qual- 
itative remarks here are useful. The essential features 
of the FEDVR method are 1) to subdivide space into a 
variable number of elements of different length and 2) to 
use a flexible number of DVR basis functions in each ele- 
ment to accommodate the local spatial variations of the 
wavefunction. The basis functions in each element are 
chosen to be the Lagrange interpolation polynomials at 
the Gauss-Legendre-Lobatto quadrature points in that 
element. The Lobatto quadrature rule, which forces a 
common quadrature point at each boundary, enables the 
construction of a "bridge" basis function which ensures 
wavefunction continuity (communication) across the el- 
ements. The "bridge" function is simply a linear com- 
bination of the last DVR function in region i and the 
first DVR function in region i + 1. Note, that the or- 
der of the DVR basis in each element need not be the 
same. In addition, the use of a DVR basis in each ele- 
ment, means that matrix elements of the potential are 



simply the value of the potential at the Gauss-Legendre- 
Lobatto quadrature points in that clement. The compo- 
nents of the kinetic energy matrix in each element are 
given by analytic expressions involving the derivatives of 
the Lagrange interpolating polynomials evaluated at the 
quadrature points. First and higher order derivatives of 
the basis set are not and need not be continuous across 
the elements for the numerical procedure to be conver- 
gent in the limit of a complete basis. The procedure 
enables us to capture both the short and long range fea- 
tures of the molecular wavefunctions with a reasonably 
small basis (~300 functions/channel). We thus have a 
large number of points in the potential region, where 
the wavefunction oscillates rapidly, and an increasingly 
sparse grid, as one gets deeper into the asymptotic re- 
gion where the behavior of the wavefunction is governed 
purely by wavelength (energy) considerations. Finally, 
an appealing attribute of the FEDVR basis for R-matrix 
calculations is that only a single FEDVR function con- 
tributes to the R-matrix elements at the boundary. This 
is a consequence of the fact that the DVR functions are 
either zero or one at the FEDVR points. 



V. RESULTS AND DISCUSSION 

In this section we present results of calculations within 
our two-channel model, which fully illustrate the reso- 
nance scattering and molecular physics of the Feshbach 
resonance. In Table [I] the parameters defining the mini- 
mal two-channel model description and the set of phys- 
ical parameters on which it is based are given for sev- 
eral Feshbach resonances of experimental interest. We 
present a comprehensive comparison of our two-channel 
results with those of a full coupled channels calcula- 
tion only for the specific example of the Feshbach res- 
onance occuring at i?o=202.15 G in fermionic 40 K in the 
|F = 9/2, nip = -7/2) and \F = 9/2,m F = -9/2) hy- 
perfine states. However, we have checked that our model 
is in agreement with more elaborate calculations for all 
the cases quoted in the table. 





Physical parameters 


Model parameters 


Species 


Ca (au) 


fflbg (au) 


N s 


Bo (G) 


AB (G) 


Afi (MHz/G) 


e (au) 


a (au) 


Wo (au) 


E th (au) 


4U K 


3897 


172.06 


27 


202.15 


7.788 


2.35 


i.63i-io~ d 


9.177 


0.221 


1.192-10"' 


4U R 


3897 


172.15 


27 


224.21 


8.12 


2.35 


1.631-HT a 


9.177 


0.226 


1. 112-10-' 


s:, Rb 


4700 


-480 


35 


155.199 


10.641 


-3.345 


1.124-10 -3 


10.075 


0.203 


3.283-10-' 



TABLE I: Physical parameters for a Feshbach resonance 
and the corresponding parameters describing the two-channel 
model. The units of the C% coefficient is Energy-Length 6 , 
with 1 au = 9.57344- 10 -26 J nm 6 . A frequency of 1 MHz 
corresponds to a temperature of 48 /iK. 



The essential molecular physics of the Feshbach reso- last two bound states of the coupled channels problem are 
nance is illustrated in Fig. [21 where the energies of the 
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plotted (solid lines) . The energies of the uncoupled dia- 
batic states corresponding to the bare Feshbach molecule 
with slope A/i and the last bound state of the open chan- 
nel at energy E_]_ are indicated by the dashed lines. In 
this case the coupling shifts the Feshbach resonance posi- 
tion downwards from the zero crossing of the bare Fesh- 
bash state at B = J3 res =211.39 G to B , where the bind- 
ing energy of the last bound state at B=0 vanishes, and 
this state enters the continuum. The last bound state of 
the open channel is seen to play a pivotal role, as it inter- 
acts strongly with the bare resonance state, leading to an 
avoided crossing. Hence the dressed molecules produced 
in an adiabatic sweep of the magnetic field across the 
Feshbach resonance, starting at B > Bq, will evolve into 
the last bound state of the open channel potential (ma- 
genta curve) for B <C B , not the bare resonance state 
in the closed channel. Instead it is the dressed state with 
energy E_ 2 (green curve) which tracks E [es (B) in this 
regime. At fields much larger than Bg the energy of this 
state approaches E^\. This agrees with earlier observa- 
tions by other authors |14| . For other resonances, such 
as the two in 6 Li, the last bound state of the open chan- 
nel lies much deeper, and has no effect on the Feshbach 
molecules. 

The open circles in Fig. [21 give the results of a full 
coupled channels calculation. We observe that our sim- 
ple two-channel model gives a faithful representation of 
the dressed states. The figure also shows the discrete R- 
matrix representation of the scattering continuum. At 
low energies they bear some semblance to the eigenstates 
in a spherical box with zero derivative on the bound- 
ary, though no such correspondence can be established 
at higher energies. Their spectral density is therefore 
tied to the size of the box. In particular, the first contin- 
uum state is not at exactly zero energy, but is displaced a 
small positive amount, representing the zero point energy 
in the effective box potential. The scattering resonance 
is visible as a scar across the continuum part of the spec- 
trum. 

We remark that while the energy of the second to last 
adiabatic bound state approaches that of the bare res- 
onance state far below the resonance, it is in general 
affected by the coupling to deeper bound states of the 
open channel potential and in fact never becomes linear 
in the magnetic field. It is therefore important to obtain 
the magnetic moment difference by fitting the position of 
the scattering resonance for fields above Bo- 
Figure [21 gives additional details of the behavior of the 
dressed states as the magnetic field is varied. The figure 
displays the weight of the dressed state wavefunctions in 
the closed channel, defined as 

Z„ {B) = J dr\ijj 2 , n (r)\ 2 . (33) 

In accordance with the bound state spectrum in Fig. [21 
the last adiabatic bound state is a pure open channel 
state in both the high and low field limit. Near the Fes- 
hbach resonance it acquires some closed channel char- 
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100 200 300 400 

BIG] 

FIG. 2: Bound states and the discrete R-matrix representa- 
tion of the scattering continuum for the 40 K Feshbach reso- 
nance at Bo = 202.15 G. The E-i (magenta line) and E-2 
(green line) bound states of the two-channel model as counted 
at B = 0, are compared with the results of a full multi- 
channel calculation (open symbols). At the resonance po- 
sition B = Bo the number of bound states changes by one. 
The low energy R-matrix eigenstates with E > (blue lines) 
correspond to the eigenstates of a spherical box with vanish- 
ing derivative at the boundary. At higher energies no such 
correspondence exists (see text). The bare resonance state 
and the last bound state of the open channel potential with 
energy E_l = 8.90 MHz are indicated by the dashed lines. 



acter, but it remains dominated by the open channel. 
In contrast, the bound state with energy E^ 2 changes 
from the open channel state with energy E^\ to the bare 
resonance state in the closed channel, as the magnetic 
field is changed from above to below Bq [3(j ■ Again, the 
results of a coupled channels calculation are shown for 
comparison, and are found to confirm the simple model. 
For the coupled channels calculation Z n represents the 
sum of contributions from all the closed channels cou- 
pled with the open channel. The variation of Z n with B 
for the lowest lying states in the continuum shows that as 
the resonance energy increases with increasing magnetic 
field, the discrete Feshbach state becomes less diluted by 
the coupling to the continuum. 

To extract information about the low energy scatter- 
ing in the two-channel model, we find the scattering 
phase shift as a function of energy using 1|29[1 . Fitting 
to the effective range expansion k cot 5(E) = —l/a(B) + 
r e (B)k 2 /2 then gives the magnetic field dependent scat- 
tering length and effective range. These are plotted in 
Fig. and the scattering length is seen to be in agree- 
ment with (Q, as expected since the two-channel model 
has been constructed to reproduce the salient features 
of this expression. For reference the effective range is 
also plotted. For most of the magnetic field interval 
shown r e remains small and positive, but it diverges at 
B = B + A5=209.93 G, where the scattering length 
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FIG. 3: Closed channel fraction for the energy curves in Fig.|3 
Note that for the molecular state created by a downward 
sweep across the Feshbach resonance the population in the 
closed channel never exceeds 8%. The closed channel popula- 
tion for the last two bound states resulting from a CC calcula- 
tion are shown for comparison. Continuum curves with peak 
positions at higher magnetic fields correspond to R-matrix 
eigenstates with a higher energy. 



vanishes (IfHsij. 
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-1000 
-1 500 
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FIG. 4: Low energy scattering parameters for the two-channel 
model for the same resonance as above. The scattering length 
(blue open circles) and effective range (dashed line) found 
by fitting the energy dependent phase shift to the effective 
range expansion. The scattering length is compared with the 
resonance expression Q, shown by the solid line. 

The energy width of the quasi-bound dressed Feshbach 
state as determined by (|14|) is plotted in Fig. [5] As the en- 
ergy approaches zero T(i?) oc k as required by the Wigner 
threshold law. Using l|13l) it follows that 

limr( J B) = 2a hg A(iABk. (34) 

This Wigner threshold regime limit is also indicated in 
Fig.0 



FIG. 5: Resonance width for 4U K as a function of the scatter- 
ing energy on a, logarithmic scale. The two-channel expression 
l|14[l (solid line) is compared to its limit in the Wigner thresh- 
old regime 134H (dashed line). 



To illustrate in greater detail how the scattering 
evolves with energy and magnetic field we plot in Fig. [B] 
the energy dependent T-matrix, T = I — S, for several 
different values of the magnetic field. We only plot the 
purely open channel component of T, which is related to 
the phase shift via 

|Tn(£;,B)| 2 = 4sm 2 (5(E,B)). (35) 

The total phase shift may be written as a sum of a back- 
ground and the resonant contribution 0, 

S(E)=5 hg (E) + 5 res (E,B). (36) 

For the resonance under consideration here we may take 
the background term to be independent of the magnetic 
field, but this is not always the case (see e.g. For 
each magnetic field the figure compares the energy vari- 
ation of both the real and imaginary parts of Tn calcu- 
lated within the two-channel model of this paper and full 
coupled channels calculations. For completeness we also 
plot for both models the resulting |Tn| 2 . At low values 
of B the Feshbach resonance has almost no influence on 
the scattering, which is dominated by the background 
contribution. Right on the Feshbach resonance the scat- 
tering near threshold is strongly affected by the coupling 
to the closed channel bound state. In the last four panels 
the energy of the bare resonance state is indicated with 
a vertical dashed black line. When it lies in the contin- 
uum it gives rise to a scattering resonance, the shape of 
which depends value of the background phase shift in the 
neighborhood of the resonance position. At B — 218 G 
the resonance is broad (corresponding to a finite value 
of Im(Tn)), and its position does not line up with the 
energy of the bare resonance state. At higher fields the 
resonance follows the bare resonance state. 

The last three panels demonstrate the interplay be- 
tween the two different contributions to the scattering by 
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plotting Tn at magnetic field values, where the resonance 
occurs at energies such that the background phase shift 
modulo ir is close to 0, 7r/4, and 7r/2, respectively. As the 
scattering energy is increased through the resonance the 
phase shift increases by approximately tt. Correspond- 
ingly, the values of the real and imaginary parts of T\i 
reflect the value of the background phase shift on reso- 
nance. This gives ReT n (B = 243 G) w 2, ReT n (B = 
275 G) w 1, and ReT n (B = 317 G) w for the real part 
of the T-matrix on resonance, while the imaginary part 
at the same magnetic fields is lmTn(B = 243 G) w 0, 
ImT u (B = 275 G) » -1, and ImTn(B = 317 G) « 0, in 
good agreement with the numerical findings. 

VI. CONCLUSIONS 

We have presented a recipe for constructing a sim- 
plified two-channel model of an isolated Feshbach res- 
onance. Given a set of five physical parameters, which 
can be obtained from either experiments or exact cal- 
culations, the parameters characterizing our model can 
be derived without ambiguity. Within the R-matrix for- 
malism, the bound state part of the spectrum as well 
as the states scattering at arbitrary energies in the con- 
tinuum are obtained from a diagonalization of the two- 
channel Hamiltonian (including the Bloch operator). The 
FEDVR is particularly well suited to represent the model 
Hamiltonian, since the combination of short range molec- 
ular potentials with long wave length scattering requires 
a spatial grid which encompasses two widely separated 
length scales. 

We have provided the parameters defining the two- 
channel model for several Feshbach resonances used in 
current experiments, and made direct comparison be- 
tween the predictions of our model with the results of ex- 
act calculations. This comparison shows that over a large 
range of energies and magnetic fields the two-channel 
model reliably reproduces the results of full coupled chan- 



nels calculations. The usefulness of the two-channel 
model lies in its simplicity. In particular, the magnetic 
field enters only as a parameter in the closed channel po- 
tential. Hence the stationary two-channel Hamiltonian 
provides an excellent starting point for studies of prop- 
agation in the presence of a time dependent magnetic 
field. Such time-depende nt p roblems are also well suited 
for the FEDVR method 

We end by contrasting the two-channel scheme of this 
paper with simplified models of Feshbach resonances 
developed by other authors. There are examples of 
a number of similar two-channel models in the litera- 
ture [HmHIllEl. The new elements of this work 
are the use of pseudopotentials, which gives as realis- 
tic a representation of the molecular potentials of the 
full coupled channels problem as possible, and a de- 
tailed prescription for generating them. Specifically, we 
choose model potentials which support the same number 
of bound states as the entrance channel in the full prob- 
lem. We do not anticipate that the additional details of 
the model presented here will lead to fundamentally dif- 
ferent results. But it does lead to a more accurate rep- 
resentation of the dressed state wavefunctions. Similar 
remarks apply when our model is compared with those in 
which even simpler pseudopotentials are introduced, such 
as a well with a barri er p3l | , two coupled wells 0, El3 or 
zero range potentials |36( . These approaches have the ad- 
vantage that they are amenable to analytical treatment, 
but they give no details of the molecular physics at short 
range. 

The virtue of a model description is measured on its 
ability to reproduce the available data effectively, while il- 
luminating the important underlying physical principles. 
We believe that our approach coupled with an R-matrix 
analysis implemented within a FEDVR scheme consti- 
tutes an efficient and transparent procedure for gener- 
ating the bound state and scattering properties of the 
coupled channels Feshbach problem. 
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FIG. 6: The energy variation on a logarithmic scale of the 
open channel component of the T-matrix for specified mag- 
netic fields. The real and imaginary parts of Tn(E) are shown 
in blue and red, respectively, while |Tn| 2 is plotted in green. 
For each magnetic field the results of a full coupled channels 
calculation (open symbols) are compared with those of the 
two-channel model (full lines). For fields higher than B Tea the 
energy of the bare resonance state is indicated by the vertical 
dashed line. The first panel corresponds to background scat- 
tering, while in the second panel the magnetic field is tuned 
to the Feshbach resonance. For B = 243 G, B = 275 G, and 
B = 317 G the background phase shift (modulo n) is close 
to 0, 7r/4, and tt/2, respectively (compare with the plot for 
B = 100 G) 0. 



